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The cross section of reaction e + e~ — > ir + ir~-K + ir~ is calculated for energies 0.65 < y/s < 1 GeV in 
the framework of the generalized hidden local symmetry model. The calculations are compared with 
the data of CMD-2 and BaBaR. It is shown that the inclusion of heavy isovector resonances p(1450) 
and /o(1700) is necessary for reconciling calculations with the data. It is found that at y/s « 1 GeV 
the contributions of above resonances are much larger, by the factor of 30, than the p(770) one, and 
are amount to a considerable fraction ~ 0.3 — 0.6 of the latter at y/s ~ m p . 
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Among chiral models aimed at the description of in- 
teractions of the pseudoscalar mesons with the low lying 
vector and axial vector ones, see the review [l[ and refer- 
ences therein, the most elegant is the generalized hidden 
local symmetry (GHLS) model @. It relates all cou- 
pling constants to only the pion decay constant / ff and 
gpnir , and accounts for anomalous processes in a way that 
does not break low energy theorems. Strikingly, but this 
very popular model was not scrutinized in the processes 
with sufficiently soft pions where one can rely on the 
tree approximation. The purpose of the present paper is 
to fill this gap by plotting the e + e~ — * n + n~ n + n~ re- 
action cross section in the GHLS model and comparing 
the results with available data CMD-2 Q and BaBaR 
[H . When so doing, we use our recent calculations of the 
p — > An decay amplitudes [H, [|| to account for the res- 
onant production e + e~ — > p —* n + n~n + n~ . Note that 
excitations curves in [5j do not include the a\n intermedi- 
ate state [6[ nor the contact non-resonant contributions 
e + e~ — * 7* — * pirn — ► An, e + e~ — > 7* — > a\n — > An 
whose explicit form is found here. 

The ingredients for the amplitude with the resonant p 
meson are given in [5, (|. The Lagrangian of the direct 
photon coupling is 
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where g — g^^, and A^, a p , 7r stand for the pho- 
ton four- vector potential, ai(1260), n meson field, re- 
spectively. Boldface characters refer to isotopic vectors. 
Given are only the terms necessary for the n + n~n + n~ 
final state, and the contributions of the second order 
in electric charge e are neglected. Note that the con- 
tact 7* — > n + n~ and 7* — > n + n~n + n~ vertices cannot 
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be simultaneously eliminated in HLS, while the contact 
7* —* n + n~ vertex is eliminated in HLS by the parameter 
choice @. 

It is suitable to represent the energy dependence of the 
e + e~ — ► 7r + 7T~7r + 7r - reaction cross section in the form 

, v _ l2nm p T pe+e -{m p )Tf^{ S ) 
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where the leptonic width of the vector meson V on the 
mass shell looks as 
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and s — q 2 is the total energy squared in the center-of- 
mass system. The function T p ^ 4v (s) in ([2|) is evaluated 
with the effective p — > An decay amplitude M?^ 47r = 
M eff + + _ _ which includes both the resonant con- 

Pq~ >7r g 1 7r q 2 ^ 

tribution e + e~ — * 7* — * p — ► n + n~n + n~ and the contact 
one e + e~ — > 7* — > n + n~n + n~ . In the lowest order in 
electromagnetic coupling constant this amplitude is given 
by the expression 
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where e p stands for the polarization four-vector of the 
virtual p meson, and A a = A a (q 4 , q 2 , q 3 , 94), a = 1, 2, 3, 4 
are dimcnsionless invariant functions. A\ = — 1 + (1 + 
P 3i )Bi, where 
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The notations are: P a b is the operator interchanging the 
pion momenta q a <-> 9b, -D pQ b = D p (q a + qt) is the inverse 
propagator of p meson with the invariant mass squared 
(9a +9b) 2 , 



Dp(q) = m\ \-q 
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see (3.3)-(3.5) in Q for T p (y/f). The terms oc D p (q) 
in ([5]) refer to the contact terms generated by |T]). (P, Q) 
stands for invariant scalar product of two four-vectors P 



and Q, D n (p) 



p is the inverse propagator of 



pion, m n and m p are the masses of charged pion and 
p(770) meson taken from A 2 is obtained from Ai 
by interchanging q\ <-> (72, A3 is obtained from Ai by 
simultaneous interchanges 91 <-> 93, 92 94 followed by 
inverting an overall sign, and A4 is obtained from A 3 by 
interchanging q 3 «-> 94. The form of the ai propagator 
■D^ with the energy dependent width is given in 6]. 
Here T ai 7^ should be taken into account because yfs = 
1 GeV is close to m ai = 1.23 GeV (a PDG value Q) 
or to m 0l = y/2m p = 1.09 GeV given by Weinberg's 
relation. We use the approximate expression for r oi (m) 
which interpolates the curve in Q in the range 3771^ < 
m < yfs — m.^ , yfs < 1 GeV. 

The resonant contribution 7* — ► p — > it + n~ n + tt~ in 
(|4"| respects the requirement of chiral symmetry in that 
it vanishes at the vanishing momentum q ap — > (a = 
1,2,3,4) of any final pion, provided 777^ = 0. However, 
the terms due to the direct 7* — * 7r + 7r~7r + 7r~ contribu- 
tion do not vanish in the above limiting cases. This is the 
consequence of the breaking of conservation of the axial 
current by electromagnetic field, dpj^ A — eApe 3a bjp[ a 
upon neglecting the term oc m 2 . One can show that the 
terms in Q surviving in the limit q ap — * 0, correspond 
to the matrix elements of the above divergence of axial 
current. 

The results of evaluation of the e + e~ — > 7r + 7r~7r + 7r~ 
reaction cross section in GHLS model are shown in Fig. 1. 
The curves are obtained in the case ?77 Ql = 1.23 GeV; the 
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1.09 GeV look qualitatively 
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FIG. 1: The energy dependence of the e e — * ir ir ty n 
reaction cross section in the generalized hidden local symme- 
try model; m ai = 1.23 GeV. The data are CMD-2 |g] and 
BaBaR [j. "HLS" refers to the case of no 01, no contact 
terms, "GHLS" does to one with both ai meson and contact 
couplings (fTJ). "GHLS, no contact terms" refers to the model 
without contact terms. 



the same. One can see that the model is unable to re- 
produce the magnitude of the cross section at energies 
yfs > 0.8 GeV. Let us include the contributions of heavier 
resonances p 1 = p(1450) and p" = p(1700) trying to ex- 
plain the cross section magnitude at yfs > 0.8 GeV, with- 
out invoking the higher derivative terms in the effective 
lagrangian. We choose the simplest parametrization con- 
sisting of the Breit-Wigner resonance shape with the con- 
stant widths and masses m p i — 1.459 GeV, T p i = 0.147 
GeV, m p „ = 1.72 GeV, I> = 0.25 GeV taken from 
and neglect the p(770) — p(1450) — p(1700) mixing due to 
their common decay modes. This approximation results 
in no qualitative difference in the role of heavy resonance 
at yfs < 1 GeV as compared to more sophisticated mod- 
els with mixing. We also adopt the assumption of a\-K 
dominance in the p 1 ', p" — ► Air decay dynamics Q, but 
modify it to include the requirements of chiral symmetry. 
Then taking into account the p' , p" resonance contribu- 
tions results in the factor 
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multiplying the right hand side of ([2]), where Dy(q) = 
my — s — imyTy, V = p',p", s — q 2 . Free parameters 

i and ; 
of x p i is that 



x p i ana x p » are found from fitting the data. The meaning 
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analogously for x p " , where g^y = erriy / fy is the photon- 
vector meson V transition amplitude, fy is related with 
the leptonic width ([3]). Since p and p' are assumed here 
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TABLE I: Table 1. The results of fitting CMD-2 data 

xViVd.o.f m ai [GeV] 



1 -27.5 ±1.5 =0 15.4/10 1.23 

2 =0 -46.2 ±2.5 15.4/10 1.23 

3 96.8 ± 1.5 -208.7 ±2.5 14.5/9 1.23 

4 -17.8 ± 1.0 = 15.7/10 1.09 

5 =0 -30.1 ±1.5 15.4/10 1.09 

6 72.5 ±1.0 -151.9 ±1.6 14.7/9 1.09 



TABLE II: Table 2. The results of fitting BaBaR data [j. 

xV^d.o.f m 01 [GeV] 



1 -25.2 ± 0.9 = 32.6/16 

2 =0 -44.0 ±2.1 29.3/16 

3 273.2 ±1.4 -514.5 ±2.3 11.2/15 

4 -15.8 ± 0.8 = 35.0/16 

5 =0 -27.7 ±1.3 31.8/16 

6 198.5 ± 1.0 -370.1 ±1.5 11.2/15 
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FIG. 2: The results of fitting the CMD-2 data 0. "GHLS" 
refers to the model without p' , p" . 
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FIG. 3: The same as in Fig. 2, but for the BaBaR data [J]. 



to have the similar coupling to the state ai7r, the ratio 
([5]) is constant. The complex function r(s) in ([7]) is the 
ratio of the amplitude with the intermediate a\ meson 
to one with no a,\ contribution. It approximately takes 
into account the a\ix dominance in the four pion decay 
of heavy isovector resonances and is precalculated for the 
CMD-2 and BaBaR Q data points yfs < 1 GeV: 
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state only, while r*_^° ai 



_ r eff,no ai ^ ig the effective 



width of the same decay including all the contribution 
mentioned above except the aiir one. The approxima- 
tion @ corresponds to the averaging over four pion phase 
space necessary to evade unacceptably long time in the 
fitting procedure. 

The results of fitting the CMD-2 data are given in Ta- 
ble 1. The curves corresponding to the fit variant 3 with 
p' and p" resonances are shown in Fig. 2. This variant 
is indistinguishable from the variants with the single p' 
(variant 1) or p" (variant 2), both resulting in the same 
curves as the dashed one shown in Fig. 2. Variants 4 — 6 



correspond to the fits with the mass m ai 



,V2= 1.09 



GeV and result in the same corresponding curves not 
shown here. The quality of fit is not quite good. Nev- 
ertheless, we quote the contribution of the sum p' ± p" 
(variant 3) or p' (variant 1) and p" (variant 2) relative 
to the case of pure GHLS contribution (dotted line in 
Fig. 2) to be 0.3 at y/s w m p and 32 at y/s — 1 GeV. 
These numbers refer to the case m ai — 1.23 GeV. The 
case m ai — 1.09 GeV results in almost the same figures 
for above ratios. 

The results of the similar analysis of the BaBaR data 
[3] are presented in Table 2. Contrary to the previous 
case, here the variants with the single additional heavy 
resonance give a bad description. The fit chooses two de- 
structively interfering p' and p" resonances each coupled 
to ai7r much strongly than in the variants of the single 
heavy resonance. The curves shown in Fig. 3 refer to 
variant 3 in Table 2 with m ai — 1.23 GeV. The contri- 
bution of the sum p' + p" ( in variant 3) or p' (variant 
1) and p" (variant 2) relative to the case of pure GHLS 
contribution (dotted line in Fig. 3)is found to be 0.6 at 
y/s « m p and 30 at yfs — 1 GeV. As in the case of the 
CMD-2 data, here the variant 6 with m ai — 1-09 GeV 
results in practically the same corresponding curves and 
ratios. 

Our conclusions differ from the result of the works 
IE B [HI all claiming small or even absent contribution of 
heavy resonances. We attribute this disagreement to the 
difference among the models used in the present analy- 
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sis and in works [!, H, H, [T(| ■ The works 0, [T(| exploit 
non-chiral invariant effective Lagrangians. The work Q 
is based on chiral amplitude with three unknown param- 
eters. No central values nor their errors are given in 
order to assess independently the quality of approach [9j . 
The effective vertex aipw used in that work refers to the 
higher derivative contribution, while there exists a lowest 
derivative one used in the present work, see Q . The con- 
tact 77r + 7r - vertex is present in the intermediate state of 
the amplitude in 0] . The apparent violation of the vec- 
tor dominance of the pion form factor could be evaded by 
adjusting arbitrary constants in in [9j only assuming the 
vanishing of the p meson width which is inappropriate in 
the energy range where the p width is essential. 

Thus, the simplest variant of GHLS model with the 
minimal number of derivatives fails to explain the cross 



section of the reaction e + e~ — > 7r + 7r~7r + 7r~ at energies 
0.8 < y/s < 1 GeV. One possible way out this diffi- 
culty by including heavy resonances p', p" is studied here. 
GHLS model is based on the nonlinear realization of chi- 
ral symmetry. It would be desirable to readdress the 
present issues in the frame work of the chiral model of 
the vector and axial vector mesons based on the linear 
cr-model. This task is necessary in order to evaluate the 
robustness of the figures characterizing the contributions 
of heavier resonances towards various model assumptions 
and to reveal the role of the intermediate states which in- 
clude the widely discussed scalar a meson. 
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